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^ , Abstract. We use analytical and numerical methods to investigate the equations for cohomo- 

^ ■ geneity one shrinking gradient Ricci solitons. We show the existence of a winding number for this 

system around the subvariety of phase space corresponding to Einstein solutions and obtain some 
. estimates for it. We prove a non-existence result for certain orbit types, analogous to that of Bohm 

in the Einstein case. We also carry out numerical investigations for selected orbit types. 

o' 
q 

(~| ' 0. Introduction 

» I 

^ - A Ricci soliton consists of a complete Riemannian metric g and a complete vector field X on a 

manifold which satisfy the equation: 

(N: (0-1) Ric(5) + ^Lx<7 + |5 = 

ijr^ . where e is a real constant and L denotes the Lie derivative. The soliton is called steady if e = 0, 
0^ \ expanding if e > 0, and shrinking if e < 0. If X is a Killing vector field, the metric is actually 

Einstein, and the soliton is called trivial: we shall exclude these in what follows. For general 
background, see [Cetcj and |Ca2] . 

As well as being natural generalisations of Einstein metrics, Ricci solitons are of interest as 
generating particularly simple solutions to the Ricci flow, for they evolve under this flow just by 
its natural symmetries, that is, by diffeomorphisms and homothetic rescalings. Ricci solitons are of 
crucial importance in analysing, via dilations, singularities of the flow |Pej. The Ricci flow is also of 
physical interest as an approximation to the renormalisation group flow for nonlinear sigma-models 
[F] . Indeed one of the earliest nontrivial Ricci solitons, the Hamilton cigar [Haj, was also discovered 
in a physical context by Witten [Wi]. 

Most known examples of Ricci solitons are of gradient type, that is, we have X = grad u for a 
smooth function u. Equation (jO.ip then becomes 



in 
o 



(0.2) Ric(c/) + Hess(n) + | g = 0. 

We shall use the abbreviation GRS for gradient Ricci solitons. In this case, the completeness of 
grad u actually follows from the completeness of g, by the recent work of Z. Zhang [Zhj . 

Noncompact complete solitons may be steady, expanding or shrinking, but nontrivial compact 
solitons must be shrinking and of gradient type [Pej . In all cases Kahler examples have been 
constructed (see [Koi] . [Calj . [ACGTj . |FIK] . |DWlj . |FuWj for example, among many other au- 
thors) . In the expanding and steady cases non-Kahler examples were constructed in |DW2j , [DW3] , 
generalising examples of Bryant |Bry| , Ivey [Ivj , Gastel and Kronz |GK| . There are also many ho- 
mogeneous non-Kahler expanders (necessarily of non-gradient type), constructed by Lauret |Laj 
and others. 

In contrast, so far the only known complete nontrivial shrinking solitons are Kahler, or else are 
rigid (i.e. obtained from a product of Einstein spaces and a Gaussian soliton [P W] ) . In this paper 
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we begin a study of shrinking solitons which are not necessarily of Kahler type. We shall focus on 
solitons of cohomogeneity one, when the existence of a symmetry group with hypersurface principal 
orbits (or a suitable bundle ansatz) reduces the equations to a system of ordinary differential 
equations. This is a natural place to start, as homogeneous compact solitons are trivial. Moreover, 
the great majority of known complete nontrivial solitons are either homogeneous (in the case of 
Lauret's noncompact expanders) or obtained by a cohomogeneity one/bundle-type construction. 
(The main exceptions we are aware of are the existence result |WZhj for Kahler-Ricci solitons on 
toric varieties, and the generalisation [PS] to Fano toric bundles over a generalised flag variety). 

In §1,2 we show that the general cohomogeneity one shrinking soliton equations yield a winding 
number, representing the winding of the flow around the subvariety of phase space corresponding 
to Einstein solutions. We obtain some estimates for the winding. We also introduce a functional 
which is monotone in a certain region of phase space. In §3 we prove a nonexistence result for 
compact cohomogeneity one solitons with certain orbit types, analogous to a result of Bohm in the 
Einstein case |B3j . In §4 we focus on solitons of multiple warped product type. In §5 we report 
on numerical investigations for some orbit types, including those where solutions were found in the 
Einstein case |B1]. The numerics indicate that the shrinking soliton equations on compact spaces 
exhibit a high degree of rigidity; in particular we discuss several spaces where Einstein metrics exist 
|B1] but where numerical searches did not produce any evidence of nontrivial solitons. 

Acknowledgements: The second author was supported by an EPSRC doctoral grant during the 
period when some of this research was carried out. He would like to thank his PhD supervisor. 
Prof. Simon Donaldson, for his encouragement. The third author is partially supported by NSERC 
Grant No. OPG0009421 

1. Cohomogeneity one Ricci soliton equations and a general winding number 

Let us consider a cohomogeneity one gradient Ricci soliton: this means that the soliton is invari- 
ant under the action of a Lie group G with generic orbit (called the principal orbit type) of real 
codimension one. We write the principal orbit as G/K, where G is a compact Lie group and is a 
closed subgroup. We refer to |DW1] for background on the soliton equations in the cohomogeneity 
one setting. 

We denote by n the dimension of G/K, so the cohomogeneity one manifold M has dimension 
n + 1. The open and dense subset Mq of the manifold consisting of the principal orbits may be 
written as / x (G/K), where / is an open interval in M and the product structure results from 
choosing a geodesic which intersects all generic orbits orthogonally. On Mq we may then write the 
soliton potential as u{t) and the metric as 



where, for each t, gt is a G-invariant metric on G/K. Letting p be an Ad/f-invariant complement 
to t in g, we know that G-invariant metrics on G/K correspond to Adi<--invariant inner products 
on p. 

The cohomogeneity one gradient Ricci soliton equations on / x [G/K) are now 




g ■= dt^ + gt 



(1.3) 



(1.2) 




(1.4) 



d{iTL) + 5^L = 



where L denotes the shape operator of the hypersurface orbits, viewed as an endomorphism and 
related to the metric by gt = 2gtL. Also rt denotes the Ricci endomorphism, and 5^ is the 
CO differential for T*(G/X)-valued 1-forms. Note that the first equation represents the components 



COHOMOGENEITY ONE SHRINKING RICCI SOLITONS: AN ANALYTIC AND NUMERICAL STUDY 3 



of the equations tangent to G/K, the second is the equation in the normal direction, and the last 
equation represents the mixed directions. 

The cohomogeneity one manifold can have up to two special orbits, corresponding to the end- 
points of the interval I. In this paper we will assume that there is at least one special orbit, placed 
at t = 0. We may then assume that / = (0, T) where T is an extended positive real number. If 
T < +00, M is compact and additional conditions must be imposed at and T to ensure that the 
metric g and potential u extend smoothly to all of M. When there is only one special orbit, the 
cohomogeneity one manifold is non-compact, and there are further conditions which ensure that g 
is complete. (See |EWj and |DW1| for more details of the smoothness conditions.) Here we just 
like to highlight one result from |DWlj (cf Prop 3.19) which has practical consequences for the 
construction of gradient Ricci solitons. 

Assume that there is a special orbit of dimension strictly smaller than n. Then if ()1.2p holds for 
a sufficiently smooth (C^) metric and potential, then (jl.4p automatically holds. Furthermore, if 
the well-known conservation law 

(1.5) u + {-ii + tr L)u -eu = C 

holds for some fixed constant C, then ()1.3p also holds. In other words, in constructing cohomogene- 
ity one gradient Ricci solitons, it is enough to produce a solution to (jl.2p satisfying the conservation 
law which is smooth up to order 3. Note that in [EW] examples of principal orbit and special orbits 
are given for which the smoothness conditions for g have to be checked to arbitrarily high orders. 

Using (jl.3p and the trace of (jl.2p . the conservation law above can be written in the form 

(1.6) S + tr(L2) - (-n + tr if + {n - 1)^ = C + eu, 

where S is the scalar curvature of the principal orbits and the constant C is the same one as in 
(II. Sp . We note that this form of the conservation law may be viewed as the zero energy condition 
for a Hamiltonian constructed out of Perelman's energy and entropy functionals. 



Two quantities appearing in (jl.6p are of interest to us. 
The first is 

(1.7) ^:=-tt + trL, 

which is the analogue of the mean curvature trL of the principal orbits for the dilaton volume 
element e~"d/Xg. By ()1.3p . ^ is strictly decreasing in the shrinking case (e < 0). Indeed we have the 
useful inequality 

(1.8) ^ ^ i < 0- 

Furthermore, when there is a special orbit with codimension k + 1 > 2 placed at t = 0, the 
smoothness conditions imply that n(0) = and trL = j + 0{t) as t — >• 0. In the complete, 
noncompact case, the two facts above show that ^ tends to +00 as t tends to and ^ tends to —00 
as t tends to +00. In the compact case, an analogous calculation to that above shows that ^ tends 
to —00 as we approach the second special orbit at t = T. Hence in the shrinking case, if the metric 
is complete, ^ must decrease from +00 to —00 regardless of whether M is compact or non-compact. 
Denote by t^, the unique zero of ^, and introduce, as in |DW2| . [DW3| . the variable 

(1.9) PF:=^-i = _i — 
as well as a new independent variable s defined by 

(1.10) #:=^^ ' ^ 



ds dt {—ii + trL) dt 
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We shall use a prime ' to denote differentiation with respect to s. It is natural to let s = — oo 
correspond to t = 0. The zero of ^ corresponds to the blow-up time of W, which occurs at a 
finite value of s. To see this, note that equation (I1.3P can be rewritten as 

C = -tr(L2) + 1, 



or in terms of W as 

(1.11) W' = W''^{tv{L 



3 / . ,^os e 



2. 

Since e < 0, we have W' > — which shows that W blows up in finite time. 

We shall call a point where W blows up (or equivalently, where the generalised mean curvature 
^ vanishes), a turning point. For the special case of Einstein trajectories, where u is constant, 
this reduces to Bohm's notion of turning point, that is, a point where the hypersurface becomes 
minimal [Bl| . Note however that a turning point is not a singularity of the metric, just a coordinate 
singularity. One might hope to construct compact shrinkers by glueing together two trajectories at 
the turning point if the right matching conditions held. 

The second interesting quantity is 

(1.12) £ = C + eu, 

which we view both as a smooth function on M and on the interval (orbit space) /. Since one 
may add an arbitrary constant to the soliton potential, ()1.5p shows that this is compensated by 
a corresponding change in the constant C appearing in the conservation law. This is one of the 
reasons for interest in the quantity £ rather than the potential. It is then natural to introduce also 
T := ii, as u is defined up to an arbitrary constant. Note that ()1.5p may now be written as 

(1.13) £ + ^£-€£ = 0. 
We now switch to derivatives with respect to s. We have 

£' = eu' = tWT. 

Differentiating T we obtain 

T' = W u 



WiC + eu 



u 
W 

where we used (jl.Sp to get the second equality above. In other words, we have the associated 
system 

£ \' f Q eW \ f £ 



along the solution trajectories as long as the independent variable s makes sense. 
Note that the eigenvalues of the matrix are 



-1 ± Vl + ^eW^ 
2 

so in the shrinking case (e < 0) they form a complex conjugate pair with negative real part, provided 
W > 



We will now examine the system (|1.14p more closely. 
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Remark 1.15. The origin in the J-") -plane (which is invariant under the flow) represents tra- 
jectories with constant soliton potential n, that is, Einstein trajectories. 

Remark 1.16. In the expanding case e is positive, and hence the quadrants i5, < and £,J->0 
are invariant under the flow (cf [DW3J ) . In [DW3J this fact was used to establish bounds on the 
flow for the case of expanding solitons of multiple warped product type. 

In the shrinking case this is no longer the case. Indeed, £' has the opposite sign to J-, while J-' 
takes the same sign as £^ on = 0. The phase plane diagram shows that no quadrant is invariant 
and indeed the flow exhibits rotational behaviour around the origin in the (iS, J-")-plane. 



In order to analyse the behaviour of the variables £ and J- in the shrinking case, it is convenient 
to take e = — 1. We observe, from the equations (jl.l4p . that 



so £^ + J-"^ is monotonic decreasing. (Note that for a non-trivial soliton, T = ii cannot be along 
the flow.) 

We also calculate 

d F' T — T' F F T 

(1.17) iLtan-i(£:/J-) = ^V-^ = -W^ + ^^^- 

If we let Q = tan^^(iS/J^), this equation becomes 

(1.18) ^ = -I^+^sin(20). 

as I 

We are interested in estimating the winding angle around the origin of the flow in the (i5, 7-")-plane 
up to the (unique) turning point (point where W blows up, cf the above discussion). 

Let us take the J^-axis as horizontal and the f-axis as vertical. We shall consider trajectories 
emanating from a point on the positive i5-axis (i.e., = ^) since at a special orbit we need u = 0, 
that is, = 0. Our phase-plane diagram shows that the flow will initially move into the quadrant 
with f , > and can never re-cross the positive i5-axis from the right. The winding angle is 
therefore nonpositive. Analogous statements hold for trajectories emanating from a point on the 
negative <S-axis, where Q = 

Observe further that the line ^ = ^ is a reflecting barrier for the flow unless W >\. Also note 
that once W exceeds \ then Q is monotonic decreasing, i.e., the flow winds clockwise around the 
origin. 

The winding angle is 

e'ds= e'ds+ e' ds, 

J-oo Jw-^{^) 

where we have arranged for s = — oo to correspond to t = (in view of (|1.10p and the initial 
conditions). The value of the first integral is between zero and — ^, from our above remarks, as 
until W reaches ^ the flow is trapped between ^ = f and 9 = j. 

For the second integral, we can change variables to rewrite it (using (jl.lip with Q := M^^tr(L^) 
and (fTTH]) ) as 

1 /■°° ^ sin(26') 

dW + / VAtt dW 



1 g + ^w^ J I w{g + \w^) 

These integrals are bounded in absolute value by 4 and 2 respectively, just using the nonegativity 
of Q. We have therefore deduced 
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Proposition 1.19. For trajectories of the associated flow in the {£ ,J^)-plane starting from either 
the positive £ or negative £ axis, the winding angle around the origin up to the turning point is 
finite, nonpositive, and hounded below by ~{6 + j). □ 

Of course one would like solutions with winding angle — vr, so as to get = at the turning point 
Sif, which is one of the conditions needed to get a symmetric solution. The above argument shows 
that one cannot achieve this by mimicking the Bohm techniques-we cannot make the winding angle 
arbitrarily large by varying some parameter. 

Remark 1.20. We can relate £, J- to the quantities used in the paper on expanding solitons |DW3j 
as follows. Let us define 

(1.21) n:=W iiL and Q := W^£. 
We recall from |DW3j the calculation: 

(1.22) u' = = WtiL -l = n-l 

^ ^ -u + trL 

which is actually valid for general orbit types. It follows that J- = . 
We can then derive easily the equations: 

(1.23) {%-!)' =(g -I- (^ - 1) + Q, 

(1.24) Q! = eW"^ {n-l) + 2{g- Uv^'^ Q, 

which we saw in jDW3j in the special case where e was positive and the cohomogeneity one metric 
was a warped product on m factors. 

The associated system (jl.l4p described earlier may be regarded as a renormalised version of the 
above pair of equations where we replace 71 — 1 and Q by their quotients by appropriate powers 
of the variable W. Note that the subset {Ti = 1, Q = 0} is invariant under the flow of the GRS 
equations (at least while the variable s is defined), and consists of Einstein trajectories. 

2. Some properties of complete shrinkers of cohomogeneity one 

In this section we describe some properties of complete shrinking gradient Ricci solitons with 
cohomogeneity one. 

It is known that the scalar curvature of a complete steady or shrinking gradient Ricci soliton 
must be non-negative. In the compact case, steady GRS are in fact Ricci flat and shrinking GRS 
must have positive scalar curvature (cf Proposition 1.13 of |Cetcj ). For the non-compact case, the 
non-negativity follows from a theorem of B. L. Chen (Corollary 2.5 of |Ch| ) on complete ancient 
solutions of the Ricci flow. (The result was previously known assuming sectional curvature bounds.) 

For cohomogeneity one gradient Ricci solitons, the ambient scalar curvature R is given by 

(2.1) R = -2 tr(L) - tr(L2) _ (trL)^ + S, 

where S denotes the scalar curvature of the principal orbits (cf §1 of [DWlj ). Now we have 

R = -S+ (tvLf - tr(L^) - 2u trL - en 
= -C- eii - - |(n + 1) 
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where we used (ll.2p and (ll.6p respectively in the first two equalities. Now R> immediately leads 
to the upper bound 

(2.2) £<-i(n+l), 

which is strict in the compact shrinking case. 

We take this opportunity to emphasize that £ can be viewed both as a function on M and a 
function on /, and we will frequently switch from one viewpoint to the other in the following. 

Remark 2.3. For the steady case, the above inequality becomes C < with equality iff ii is 
identically zero, i.e., the soliton is trivial. 

Remark 2.4. In the case of a complete non-trivial expanding GRS of cohomogeneity one, Corollary 
2.3 of [Chj . applied to the Ricci flow generated by the expander, implies that R > — |(n-|- 1). (Note 
that our ambient manifold has dimension n + 1.) Since e is now positive, the trace of (|U.2|) yields 
AS < 0, where A is the negative Laplacian of the cohomogeneity one metric. As £^ is a function 
of t only, AS = S + (tTL)£. Using the above facts in (I1.13P we obtain £ < 0. In particular, if we 
assume without loss of generality that 'u(O) = 0, then the constant C in the conservation law must 
be negative. 

Going back to the shrinking case, recall that the smoothness conditions at a special orbit placed 
at t = imply that 

(2.5) u = n(0) + ^ u{0) + ■■■ 
and 

(2.6) trL = ^ + • • • 

for sufficiently small values of t. (Here the mean curvature trL of the principal orbits is taken 
with respect to the "outward" normal ^.) Hence, in the case of a complete shrinking GRS of 
cohomogeneity one, the conservation law (jl.Sp and (|2.2p give 

(2.7) {k + l)u(O) = C + en(0) = £{t = 0) < -^(n + 1). 
So for the initial value ii(0) we have the upper bound 

(2.8) .(0) < -5 

which must be strict at both endpoints in the case of compact shrinkers. This means that the 
trajectories of the associated flow in the (iS, J-")-plane cannot start arbitrarily high on the positive 
i5-axis if they correspond to complete solutions, compact or not. By contrast, note that there is no 
such restriction on u{Q) from the initial value problem (cf |Buj ) around special orbits. 

Note also that the upper bound £{t = Q) < — |(n-|- 1) is attained by the conical Gaussian soliton 

on M"+i. Indeed, in this case £ = -^{n + I) - (see Remark [412]) . 

Proposition 2.9. Let (M, g) he a non-trivial compact GRS of cohomogeneity one with G-invariant 
potential function u. Let [0,T] denote the closure of the interval L and £ = C + eu, where the 
constant G is that appearing in <\1.6^ . Then: 

a. £ must change sign and is a M or se-Bott function on M. 
h. £ < -f(n + l). 

c. £ has at most 4 critical points in (0,T). //to S (0, T) is a critical point then it is either a 
local maximum if £{to) > or a local minimum if£{to) < 0. 
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Proof. We note first that by ()1.13p £ can be viewed as a solution to the linear equation 

On (0,T) the coefficient ^ is finite- valued and at least C^. Since the soliton is non-trivial, £ cannot 
vanish at an interior critical point tQ. Also, <f (to) = ^^ito) holds at such a point, and the only 
non-vanishing component of the Hessian of £ (as a function on M) is that in the ^ direction. So 
the critical submanifold of £ at to, a principal orbit, is non-degenerate. This equation also implies 
the second part of c. 

Next we consider the situation at t = 0. (The situation at t = T is completely analogous.) By 
smoothness, we have £{0) = en(0) = 0. If £{0) = 0, our soliton would be trivial, as the origin 
{£, J-) = (0, 0) is invariant under the associated flow (|1.14p and represents the Einstein trajectories. 
Now the computation (]2.7p shows that n(0) = -£{0) has the same sign as £{0). In particular, 
the two singular orbits are also non-degenerate critical submanifolds of u, so that u is Morse-Bott. 
Since M is compact, £ can only have a finite number of critical points, considered as a function on 
/. Because 6 + j< |7r, we may apply Proposition II . 191 separately to each of the intervals [0, U] and 
[t*, T] to estimate the total number of critical points in (0, T). (Here as usual denotes the unique 
turning point). Indeed, when £{0) and £{T) have the same sign, the number of critical points is 
odd, and so is bounded above by 3. When £{0) and £{T) have opposite signs, then number of 
critical points is even, and so is bounded above by 4. 

The inequality in part b is just ()2.2p for the compact case. 

It remains to show that £ must change sign. The (negative) Laplacian of £ is 

A£ = £ + {tvL)£ = e£ + ^ £\ 

where we have used p.l3p in the second equality above. Hence the integral of £ is negative if the 
soliton is non-trivial. 

Suppose that £ <0 on [0, T]. By what we already proved, £ must then be negative at t = and 
t = T, and hence £ must be positive at these endpoints. Therefore, £ must have interior critical 
points. At each such point £ is negative, so that £ is positive. We obtain a contradiction at the 
first interior critical point, which must be a local maximum because £ is convex for a little while 
after t = 0. □ 

Remark 2.10. At present the only known examples of nontrivial compact shrinking Ricci solitons 
are Kiihler. These are generalizations of the Koiso examples (cf |ACGT) . |DWlj and references 
therein), and £ is monotone (no interior critical points) in these examples as a result of the Kahler 
condition. In particular, £^ is a perfect Morse-Bott function. 

In order to state the analogous result in the complete non-compact case, we need to recall 
a general result of Cao and Zhou ( |CZ) . Theorem 1.1) concerning the behaviour of the potential 
function. Adapted to the cohomogeneity one situation and with our notation, their theorem implies 
that, for sufficiently large t, 

(2.11) - |(n + 1) + ^{tV^ + C2f < £{t) = C + eu{t) < -1 (n + 1) + \{tV^e - cif 

where ci , C2 > are constants which depend on dim M = n + 1 and the geometry of the unit 
ball in M centred at a point on the singular orbit, e.g., the initial point of the solution trajectory 
associated to the soliton. 

Proposition 2.12. Let {M,g) be a non-trivial complete, non-compact GRS of cohomogeneity one 
with G-invariant potential function u. Let £ = C + eu, where the constant C is the one appearing 
in () j.gp . Then: 

a. £ must change sign and is a Morse-Bott function on M . 
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b. £ satisfies i\2.13 above. In particular, E < — |(n+l) and eventually decreases monotonically 
to —oo. 

c. £ has at most 5 critical points in (0, +oo). IftoG (0, +oo) is such a critical point then it is 
either a local maximum if £[tQ) > or a local minimum if£{to) < 0. 

Proof. The local properties of £ follow as in the compact case. The inequality (12. lip implies that 
£ tends to — oo. Suppose £ < everywhere. Since £{0) < and £ is convex for a while after 
t = 0, there must be a critical point at some to > 0. At the first such critical point we have a 
contradiction as in the compact case. So £ changes sign. 

Next choose ti such that £ < for all t > ti. So any critical point t2 > ti must have £{t2) > 
and hence must be a minimum. Since £ tends to — c«, we now see that there are actually no critical 
points in [ti,+oo), and so the number of critical points is finite. It also follows that £ < for 
t > ti. 

To get a more precise bound on the number of interior critical points, we will modify the argu- 
ments in §1 leading up to Proposition II . 191 (Note that here we do not need the work of Cao-Zhou.) 
Recall that ^ has a unique zero at t = t*. For the portion of the trajectory with t < t*, the proof 
of Proposition 11.191 gives a lower bound of — (6 + ^) for the winding angle. For the portion with 
t >t^:, we can introduce the independent variable 

(2.13) ^ ~ ^* ~ It ^^^"^ ~ 

Differentiating £,J^ with respect to a (with e = —1), we now have the following analogue of 

rift 1 

l^ = W-=- sm{2ft) 
da 2 

Note that the angle 9 is well-defined and smooth as t passes through and as t — > +oo, so does a. 
We need to estimate the integral 

de , 

— da. 

's, da 

Recall that W is now negative and increases strictly to as o" — )• +oo and satisfies the equation 

Let 0"! be the value such that W{ai) = — ^. As in §1, the integral on [s*,(Ti] is bounded below 
by —6. On [cji, +co) the net change in ft is bounded below by — tt because whatever 9{ai) equals, 
one encounters a reflecting barrier of the flow (in the clockwise direction) at or before 0{ai) — tt. 
Hence the integral over the entire trajectory is bounded below by —(6 + f ) — (6 + vr) > — (5.1)7r. 
This gives the desired bound for the number of critical points. □ 

Remark 2.14. So far the only examples of complete non-compact shrinking Ricci solitons are the 
Kahler examples of [FIK] and their generalizations (cf [DWl] ). and (up to covers) the Gaussian 
soliton and its product with a finite number of Einstein spaces. In all these cases, £ is monotone 
decreasing. Critical points of the soliton potential for complete non-compact shrinkers have been 
studied in [Na] (cf Corollary 2.1) assuming Ricci bounds. 

Next we introduce a function which has the Lyapunov property for the cohomogeneity one GRS 
equations in a region of phase space. This generalises an analogous function due to C. Bohm |B2j 
for the Einstein case. Let 

2 



(2.15) 3-:=?;n (5 + tr((L(0)) 
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where v{t) = \/ det gt , S is the scalar curvature of the principal orbit, and L^^^ := L — is the 
traceless part of the shape operator. Modulo the conservation law (II. 6p . 3" can also be written as 



2 / 11 — 1 

(2.16) J = (trL)^ + u(u-2 trL) +£ - -(n-V 

\ n 2 \ 

Proposition 2.17. Along the trajectory of a cohomogeneity one GRS of class we have 

3'=-2v^ tr((L(0))2) (^e-itrL^ 
Hence 3" is non- decreasing wherever ^ < ^trL. 

Proof. If we differentiate (|2.16p . use v = (trL)t> and ()1.3p . we obtain after some simplification 

j = 1 t;f(trL) f i(trL)2+tt2-2M(trL)+£:--(n-r 
n \ n 2 ' 

+ (trL) e u + 2tr(L2) — 

\ n n \trL n 



j = V- {ir L) [ ^— ^{iiLf - -{u-i? + 2u{iTL) -S) 



Splitting up tr(L^) as tr((L('^))^) + -(trL)^ and simplifying gives 

ytrL n y n J 

= .^(trL)(2tr((LW)2)(A_i + i 

= -2^;itr((LW)2)(e-^(trL)), 

as required. Note that we used the conservation law (II. 5p in the penultimate equality. □ 

Remark 2.18. The regularity assumption is natural since for a GRS the conservation law ()1.5p 
is usually derived from the equation for the Laplacian of du, cf Remark 3.8 in |DWlj . Conversely, 
as discussed in §1, in constructing a complete GRS, this hypothesis allows us to work with (jl.2p 
and (jl.Sp as long as there is a singular orbit of dimension strictly smaller than that of a principal 
orbit. 

3. Non-existence results 



Fundamental results of Perelman |Pe] imply that on a closed manifold all Ricci solitons are of 
gradient type and all non-trivial solitons are shrinkers. If the soliton metric is invariant under a 
compact group of isometrics, then by averaging over the group one obtains a new soliton potential 
that is an invariant function. In particular the metric of a compact homogeneous Ricci soliton 
must be Einstein. Starting from dimensions > 12, examples of compact (simply connected) homo- 
geneous spaces without homogeneous Einstein metrics are known, cf [WZ] and [BKj . Hence these 
homogeneous spaces do not admit any homogeneous soliton structures, trivial or not. 

We consider next the compact cohomogeneity one case. In this case, the orbit space is either a 
circle or a closed interval. In the former case, the fundamental group of the manifold is infinite, so 
by the work of |De] . |FGj no cohomogeneity one Ricci soliton structures exist. Let us then assume 
that the orbit space is a closed interval and the singular orbits have dimension strictly smaller 
than that of the principal orbits. Bohm showed in [B3j that for certain principal orbit types G/K 
there are no compact cohomogeneity one Einstein metrics. His hypotheses are that there exists a 
G- invariant distribution on the principal orbits G/K such that: 
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(i) the tracefree Ricci of any G-hivariant metric on G/K is negative definite on the distribution, 

(ii) the distribution is irreducible, 

(iii) the distribution has trivial intersection with the vertical (collapsing) space of the two special 
orbits. 

Many examples of such orbit types were given in [B3], §7. 

It turns out that Bohm's argument also works in the Ricci soliton case, as the extra Hessian 
terms iiL and u in (jl.2p , (jl.Sp do not play an essential role in the proof. For the convenience of the 
reader, we give a self-contained proof for the soliton case below. As a consequence, one concludes 
that the examples in |B3] also do not admit any non-trivial G-invariant Ricci soliton structures. 

Let M be a closed manifold with a cohomogeneity one smooth action by a compact Lie group 
G such that the orbit space is a finite closed interval and no singular orbit is exceptional. Let 
{g, u) be a gradient Ricci soliton structure where both the metric and potential are G-invariant. 
We may choose a unit speed geodesic 7 that intersects all principal orbits orthogonally. Then there 
are closed subgroups K C Hj C G, j = 1,2, such that G/K is the principal orbit type along the 
interior of 7 and G/Hj are the singular orbits at the endpoints of 7. Assume that the domain of 7 
is [0, r], which may be identified with the orbit space of the cohomogeneity one G-action. 

Let us choose an Adx-invariant decomposition q = t © p. Note that in this section we are not 
making any assumptions on the multiplicities of the irreducible X-summands in p. For each of the 
singular orbits G/Hj, j = 1, 2, we have a decomposition p = 5j © q^-, where t)j 5j © Sj is Ad;^- 
invariant, and is Ad/^^. -invariant. In order to address smoothness issues at G/Hi it is convenient 
to fix a G-invariant background metric f3 on G/K such that on Si it induces the constant curvature 
1 metric on the sphere Hi/K and on q^^ it induces a fixed G-invariant metric on G/Hi. (Such a 
background metric will not in general come from a bi-invariant metric on G.) Now on the open 
subset consisting of principal orbits, one can write the soliton metric g in the form dt"^ + gt where 
gt is a one-parameter family of invariant metrics on the principal orbit G/K, regarded as Adx- 
invariant endomorphisms of p which are symmetric with respect to the fixed background metric /3. 
On the other hand, the Ricci endomorphisms rt and shape operators Lt are symmetric with respect 
to gt but not necessarily with respect to /?. 

We next decompose p as a sum 

(3.1) p = Pi©---©Pr 

where each pj is itself a sum of equivalent Adi<--irreducible summands, and, for all i 7^ j, no summand 
of pj is equivalent to a summand of p^. Such a decomposition is unique up to permutation, and 
the summands pj are orthogonal with respect to /3 as well as gt- Note that any Adi^-invariant 
5(t-symmetric endomorphism of p must map pj into pj and we will denote by trj its trace on the 
summand pj. Its full trace on p will be denoted by tr as usual. 

Recall that the relative volume v{t) = \/ det gt of the principal orbits satisfies the equation 
V = {tiL)v. Consider the conformally related metrics 

_2_ 

g = v "5. 

One easily computes that 

(3.2) ~g = 25LW 
and 

(3.3) 'ti = tr^t + 2~gm, 

where L*^*^-* = L — (■^) I is the trace-free part of L. 
Consider the quantities 

(3.4) Fj := i tr.if). 
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Using (I3.2p and (13. 3p we obtain 

(3.5) Fi = tr,(5^) = 2tr,(52L(0)), 

(3.6) Fi = iTi{f) + inim-^l) + 2 iT^{~g'm). 



To compute further, we need to break the terms in Eq. (|1.2p into their trace and trace-free parts. 
Recall that the generalised mean curvature ^ is given by —u + trL. Then one obtains, using ()1.2p . 

\n 2 n \ J I 

where r'^'^^ is the trace- free part of r and S* is its trace, i.e., the scalar curvature. Note that the 
trace of Eq. p.2p yields the relation 

T) f 

S - tr(L) - e trL + y = 0. 

Hence we obtain 

(3.7) m + - r(0) = 0. 
Substituting this into (|3.6p and using (j3.5p we finally get 

(3.8) + iPi = tiiif) + tnitr^tg) + 2 tri(52rW), 

which is an analogue of the formula in Proposition 3.2 in [B3] . 
We are now ready to prove 

Theorem 3.9. Let M be a closed cohomogeneity one G-manifold as described above. Assume that 
some summand pj^ in (jg.ip is actually Adx -irreducible and that for any G -invariant metric on 
G/K, the restriction to pj^ of its traceless Ricci tensor is always negative definite. Assume further 
that Pi^ n 5j = {0} for j = 1,2. _ 
Then there cannot be any G -invariant gradient Ricci soliton structure on M . 

Proof, (cf |B3j ) Consider the function Fi^ given by ()3.4p . We claim that as t approaches 0, Fi^it) 
approaches +oo. This follows since v tends to zero (the singular orbit G/Hi has strictly smaller 
dimension than the principal orbits), and the irreducibility of p^^^ together with the fact that pj^ n 
Si = {0} imply that the endomorphism g\pi^j tends to a multiple of the identity. Similarly, we see 
that (t) tends to +oo as t tends to r. (Note that Fi^ depends on the background metric, but the 
functions defined using different backgrounds differ only by positive constants, as pj^ is irreducible.) 

It follows that Fifj has a global minimum at some point t^, G (0,r). By (13. 5p . ^(t*) restricts 
to zero on pj^^, and using this in (13. 8[) with the negative definiteness of r^*^^ on pj^^, we arrive at a 
contradiction. □ 

Remark 3.10. Examples of cohomogeneity one manifolds which satisfy the hypotheses of Theorem 
13.91 include S''^^ x {G' /K') x M3 x • • • x Mr where M3, • • • ,Mr are arbitrary compact isotropy 
irreducible homogeneous spaces, G' /K' = SU(£ + m)/(SO(£)U(l)U(m)) (a bundle over a complex 
Grassmannian with a symmetric space as fibre), and £ > 32, m = 1,2, /c = 1, . . . , [£/3] (cf [B3]). 
The significance of the spaces G' /K' is that they do not admit any G'-invariant Einstein metrics 
(cf |WZ] ) . More information about the structure of such G' /K' can be found in |B4] . cf Theorem 
B in particular. 
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4. Multiply warped principal orbit types 

In the final section §5 we are going to conduct numerical investigations of the soliton equations for 
certain relatively simple orbit types, which in the Einstein case were studied by Bohm [Blj. In the 
current section we shall therefore see how the soliton equations simplify under certain assumptions. 
We shall in particular make some remarks about the case of multiple warped products, which have 
some special features (some of the Bohm examples are warped products with two factors). 

Let us first assume that the isotropy representation is multiplicity free, that is, p decomposes 
into pairwise inequivalent real -fC- modules p]^ © . . . © p^. We let dj denote the real dimension of pj 
and let n = di + . . . + dm- Note that the summands pj have a different meaning than that in §3. 

Under the multiplicity free assumption, both the metric and shape operator diagonalise with 
respect to the decomposition pi©. . .©pm- We write the components as gf and Li = ^ respectively. 
Furthermore, the Ricci term rt also diagonalises with components rj. 

In addition to the variable W defined by (jl.Op . we introduce variables 



{-ii + tiL) Qi 
^/di 1 



(4.1) Xi 

^^'^^ ~ g^ (-li + trL)' 

Note that the non-negative quantity Q = VF^tr(L^) at the end of §1 becomes 'Y^Jj^iXf and % 

WtiL becomes YllLi Vdi^i- 

The gradient Ricci soliton system (jl.2j) - (jl.3p becomes the following equations 




(4.3) X'i = Xi\y^X^-l \ +^/dir^W^ + -{^/di-X^)w\ 

(4.4) 



(4.5) W = W l^Xf --W^ 

\j=i ^ 

Note that ()1.4p is automatically satisfied under the multiplicity free hypothesis (cf Remark 2.21 in 
[DWl] l. 

The term involving in ()4.3p can be expressed as rational functions of the Yj. Without loss of 
generality we can take W to be positive (before the blow-up time). 
In the above variables, the conservation law (II. 6p becomes: 



(4.6) C + (^^^^^ eW^ = {C + eu)W^ {= £W^) 

where C is a constant and 

^ tT(L^) tr(rt) 

^ (-n + trL)2 ^ (-n + trL)2 

In the steady case (e = 0), we recall that C is a Lyapunov function for the flow, but this fails for e 
nonzero. The first term in C is just Q = YliLi ^f- 



Let us now focus on the multiple warped product situation (see |DW2] for the steady and |DW3| 
for the expanding case). Here the hypersurface is a product of Einstein manifolds with positive 
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Einstein constants Aj. Now rj = ^ so the term involving rj in the equation for Xi is just (up 

to a positive multipHcative constant). 

Moreover, we have ti{rt) = S > and it fohows from (j4.ip and (j4.2p that 

m 
i=l 

so £ is bounded below. This can be used to give more precise information about the flow, especially 
in the neighbourhood of a turning point. By contrast, for more general orbit types, the tr(rt) term 

and the functional C involve -tpt- terms. Moreover C is indefinite. 

In the multiple warped product case we can sharpen the analysis of the functional 3" that was 
introduced in §2. As we are dealing with shrinking solitons, e is negative. 

Proposition 4.8. Let ^(t) denote a solution trajectory of the GRS equations in the multiple warped 
product case. Assume that the soliton is complete. Then the quantity rj := ^ — -trL eventually 
becomes and stays negative. Hence 9" eventually is non- decreasing, and is then strictly increasing 
unless the shape operator L becomes a multiple of the identity or we reach a second singular orbit. 

Proof. We shall assume that our soliton is non-compact; otherwise, ii is close to zero and trL is 
close to — oo near t = T, so that all the assertions hold by Proposition 12.171 

Let us consider the region {r/ < 0}. We claim that it is invariant under the forward flow of the 
GRS equations. By abuse of notation suppose ^{t) is a trajectory of the flow starting at to with 
77(7(^0))) < 0. Let ti > to be the first time when 7 leaves the region {r/ < 0}. Then 

({t^) = -{trL)\t=t,. 
n 

On the other hand 

fj = e--tr(i) 
n 

S 1 

= -tr(L^) h-^trL, 

n n 

where we have used ()1.3p and the trace of ()1.2p . At t = ti we have 

V{t,) = -i5(ti)+(^^^-tr(L2))(tO 

< -^(s+'^{tTL)A{ti)<0 
n \ n J 



since S{ti) > 0. This means that f]{'y{t)) must be positive slightly to the left of t = ti, a contradic- 
tion to the choice of ti . 

It remains to see that 'r]{'y{t)) is negative somewhere. The previous paragraph implies that once 
77(7(^1)) is zero for some ti, it must be negative to the right of ti. So suppose that r]{'y{t)) > for 
all t. Then by (fO]) . 

V 1 / ^\ ^ e 
— = -(trL) <^< -t + a 
nv n 2 

for some positive constant a for large values of t. It follows that the cohomogeneity one manifold 
has finite volume, contradicting a theorem of H.-D. Cao and X.-P. Zhu ( |Ca3j . Theorem 3.1). □ 



Remark 4.9. In the above proof the only time we used the hypothesis that the hypersurface is of 
multiple warped product type is when we asserted S{ti) > 0. 
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Proposition 4.10. Let the hypersurface (principal) orbit be a product (Mi, hi) x ■ ■ ■ x (Mm, h^) of 
Einstein manifolds with dimension di and Ric(/ii) = Xihi, Aj > 0. Assume further that Mi = S^^^ 
and hi is the constant curvature 1 metric. Then the Lyapunov function 5" satisfies 



j> nJJ a: 



di/n 



i=l 



Equality holds for a trajectory of the GRS equations iff gi/^fXi is independent of i {1 < i < m) and 
the traceless shape operator L^^^ is 0. 

Proof. By (|2.15p . we have 

/ m 2d \ ™ J \ 



,=i .91 



m 
i=l 

where we used the inequahty (cf [HLPj . p. 17) 

m _Pi 



m Pi ^ 



i=l 



for Oi, Pi > 0. EquaHty holds in this inequality iff all the Oi are equal, which translates into our 
condition since we took aj = Xj/g'j and pj = dj. 



□ 



Let us consider the 3-dimensional subset of phase space 

D := {{g,L,u) : L^^^ = 0, gil\f\i independent of i\. 



By Proposition 14.101 D is the set on which the function 3" attains its (global) minimum value. 
It is invariant under the flow of the GRS equations. In fact, if we parametrize this set by x = 

g^^J'^^,y = x/x, and z = ii, then the induced flow on D is precisely the GRS equations for the 

special case m = l,di = n, which was analysed by Kotschwar |Kot j . Recall that Kotschwar showed 
that only two trajectories in T> represented complete smooth solitons with at least one singular orbit 
(actually a point). However, when m > 1, no trajectory lying in T) represents a smooth soliton. 
When m > 1, the non-compact solution acquires a conical singularity and will be referred to as 
the conical Gaussian soliton. The compact solution has two conical singularities and will be called 
the spherical Einstein cone. On the other hand, there is a smooth non-compact complete solution 
which is the product of the shrinking Gaussian soliton on M'^^"'"^ with the remaining Einstein factors. 
(This does not represent a trajectory in D). We will call this solution the smooth Gaussian soliton, 
which is rigid in the terminology of Petersen- Wylie |PW] . For comparison purposes we list below 
some information about these special solutions. 



Example 4.11. {smooth Gaussian on 



pdi+l 



X M2 X ■ ■ ■ X Mm) The metric is given by 



dt' + t'hi+(^\ /,, + ...+ ^2A. 
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and u{t) = —\t . We further have 

Example 4.12. {conical Gaussian) The metric is given by 

\n — \ n — 1 

and u{t) = — f i^- A conical singularity occurs at t = except when m = 1. The metric is Ricci-flat 
and complete at infinity. We further have 

trL = -, ^ = -t + -, £ = t^--(n + l). 

Example 4.13. (spherical Einstein cone) The metric, which is Einstein, is given by 

dt"^ + ( T (^1^1 + • • • + ^mhm) 

\a^[n — 1) J 



where a = y^—e/2n. We also have ^ = trL = an cot(af). When m > 1 there are conical singularities 
at t = and t = n/a. 



Let us now further specialise to the case m = 2. This is of particular importance because it 
includes some of the examples for which Bohm ^Blj was able to produce infinite families of Einstein 
metrics. 

We have found, in §1, a general winding number for the flow, representing winding round the 
submanifold of Einstein trajectories. We shall now, in the case of warped products on 2 factors, 
investigate another winding number, similar to the one used by Bohm in the Einstein case (cf |Blj ). 
This counts winding around the subset D of phase space, which we introduced above. We saw that 
trajectories in D represent solutions which are equivalent to warped products on one factor. 

For convenience, we shall take the Einstein constants Aj to be — 1. We recall that D is defined 

by 




Note in particular that the intersection of 2) with the region {Ti = 1, Q = 0} (cf. Remark ll.20p 
consists (with e normalised to — 2n) of the curve parametrised by 

Vdi ^/di\/n — 1 1 

Xi = , Yi = sec t, W = — tan t. 

n n n 

This is the spherical Einstein cone solution in ()4.13p . with a = 1. 

The flow- invariant subvariety 2, of phase space given by 7^ = 1, Q = 0, = is transverse to 
D. More precisely, this subvariety meets D only at the points P± given by 

X. = }/di. . Y- = ± ~ ^ 

* n * n 

Note that P+ is the initial point (i.e. the value at t = 0) of the spherical cone and Gaussian soliton 
solutions. In fact P+ is a critical point of our equations (recall t = corresponds to s = — oo). 

Let us now linearise our equations around P+ in the invariant subvariety 2,. This means we only 
consider vectors tangent to Z, so we have 
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(the final equation follows from the condition dQ = 0). Parametrising the tangent space by xi,yi, 
we obtain the system 

/ l-n 



Xl 

yi 

The eigenvalues are roots of the quadratic 




,2 n-1 2(n-l) 



n n? 



which has discriminant (" ^)(" ^) _ \Ye deduce 

Proposition 4.14. Pj^ is a focus of the flow on the invariant variety Z transverse to D if and only 
if2<n<8. 



So this winding number behaves in just the same way as its analogue in the Einstein case in [Blj . 
where it was used to produce symmetric metrics (cf. remarks after equation (ll.lip ). However, to 
produce soliton analogues of these Bohm metrics one would need to control ii as well, probably by 
means of the winding number of §1, and as we have seen this has quite different behaviour. In the 
next section we shall investigate these and other examples of Bohm numerically in the soliton case. 

5. Numerical investigations 

The non-existence results of §3 rule out certain special principal orbit types but leave the existence 
problem open for most other orbit types. In particular they do not apply to the examples considered 
by Bohm in [Bl] , who produced infinite families of cohomogeneity one Einstein metrics on certain 
manifolds with low dimensions (5 < n + 1 < 9). 

In this section we report on some attempts to investigate the existence of cohomogeneity one 
shrinking Ricci solitons using numerical methods, focusing on the Bohm spaces and some related 
examples. The numerical methods employed are relatively simple, but do show some interesting 
results. 

We first consider an example that is known to admit a non-trivial (i.e. non-Einstein) shrinking 
Ricci soliton, namely CP^jjCP . The soliton in this case is a f7(2)-invariant Kahler metric that 
was constructed independently by Koiso |Koij and Cao [Calj . This metric provides an important 
test case for any "soliton-hunting" programs, and we shall henceforth refer to it as the Koiso-Cao 
soliton. 

For the more complicated spaces considered by Bohm we are able to replicate his Einstein metrics 
from the numerical search, but do not find any examples admitting non-trivial shrinking solitons. 
We also examine some examples of Bohm type but with dimension above the range where his 
Einstein existence results apply. In this case too we do not find any solitons. The results seem 
to hint that the cohomogeneity one soliton equations in the shrinking case exhibit a high degee of 
rigidity and instability. 

Let us now turn to the analysis of the equations. We are interested in solutions to equation 
(|0.2p . Following Bohm |B1] . we consider situations where the space of invariant metrics on the 
hypersurface is 2-dimensional. The cohomogeneity one metric therefore involves two functions, 
and the soliton equations, being a system of second-order equations in the metric and the soliton 
potential, are a dynamical system in M^. As remarked in §1, the conservation law (jl.Sp may be 
viewed as a constraint and plays a vital role. 

More precisely, we can consider a principal orbit G/K and take one of the special orbits to be 
G/H, where 

with pi,p2 being inequivalent irreducible K-modules and f) = t © p]^. As in §4, we denote by di 
the dimension of pj. To ensure smoothness at the special orbit we need H/K to be a sphere S*^^. 
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Note also that Q := G/H is isotropy irreducible and hence Einstein; we shall denote its Einstein 
constant by Cq > 0. Our cohomogeneity one metric is now given, with respect to an invariant 
background inner product i? on © by 

dt^ +gt = df + f{tfB\^^ + h{tfB\p^ . 

Note that the fibration H/K — )• G/K — )• G/H now becomes a Riemannian submersion with respect 
to gt on G/K and the metrics given by f{tfB\p^ and h{tfB\)p^ on H/K and G/H respectively. 

Examples of this situation include the case when G/K is a product of two isotropy irreducible 
spaces with one being a sphere. Another case, familiar in the Einstein situation from the work of 
Berard Bergery [BBJ, is when di = 1, so G/K is a circle bundle over G/H. 

As usual, the same equations arise in certain more general situations which are not strictly of 
cohomogeneity one type. For example in the Berard Bergery situation, the hypersurface can be a 
circle bundle over an arbitrary Fano Kahler-Einstein base, not necessarily homogeneous. 

The Ricci soliton equations (II. 3p and (II. 2p specialize to: 

(5.1) -di^-d2- + u+^ = 0, 

where the quantity is the norm-squared of the O'Neill tensor for the Riemanian submersion 

H/K — )• G/K G/H mentioned above. It is a constant that depends only on the topology of the 
bundle. When G/K is a product of the two factors H/K = S^'i and G/H then ||»4|p is zero. We 
recall the conservation law (jl.Sp . given in this setting by 



(5.4) u + u ^diy + d2^j -tt^-eu = C. 

Note that by changing the value of u{0) we may fix the constant C to be 0, and by a homothety of 
the metric we may alter the value of e. 

To carry out our numerical study we shall consider the system ()5.2p -( l5^ and introduce new 
variables {zi, Z2, zs, Z/i, z^, zq) := (/, /, h, h, n, ii). We then obtain the system 



(5.5) 


Zl 


= Z2, 








(5.6) 


Z2 


= -{dl 


zl ^^Z2Z,^dl-l^ 
Zl Z'i Zl 




1^^ + 22^6 + o^l 


(5.7) 




= Z4, 








(5.8) 


Zi 


= -{d2 


l/^ 2 II. 

ZZ Zl Z-i 




e 

+ 24^6 + -Z3, 


(5.9) 


Zb 


= Zq, 








(5.10) 


zg 


= -ZQ ( 


' , ^2 , , , 2 , 

dl 1- d2 — ] + Zq+ ez5. 

^ Zl Z^J 







The algorithm used to solve the above system numerically is the Runge-Kutta method. A good 
account of it can be found in [Butj , and we include a brief discussion here for completeness. Given an 
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ODE y{t) = F{t,y) with initial condition y(0) = yo, the Runge-Kutta sequence of approximations 
is given by: 

yn+i = ?/n + ^(ki + 2k2 + 2k3 + k4), 

D 

tn+1 = + h, 

where 

ki = F{tn,yn), 

1 1 

k2 = F(t„ + -h,y„ + -hki), 

kg = F(t, + ih,^hk2), 
k4 = F(t„ + h,hk3). 

The error at each step is 0(|h|^) which gives an accumulated error of 0(|h|^). 

The smoothness conditions at the singular orbit require that (zi, Z2, zs, z^, Z5, zq) = (0, 1, h, 0, u, 0) 
where h is an arbitrary positive constant and u is an arbitrary constant > — (n + l)/2 (by (j2.2p ). 
Note that this means that the initial Runge-Kutta step is not defined. We get around this problem 
using the standard trick of constructing a power-series solution about a small neighbourhood and 
then running the algorithm from this solution. 



5.1. Numerical results for CP^jiCP . The metrics we consider here are C/(2)-invariant and the 
principal orbits are Berger spheres = U{2)/U{1), so we have di = 1 and Q = S'^ with Cq = 1. 
The constant = 1. As remarked in the introduction, there is a non-trivial Kahler-Ricci soliton 
on this manifold due to Koiso and Cao. There is also a Hermitian, non-Kahler, Einstein metric 
due to Page [Pa]. The Page metric is constructed explicitly in our variables in [Kod] , and in order 
to check our algorithm we use the same normalisation and so take e = —7.46562. 

In order to define a smooth metric on the blowup of CP^, we need the fibre of the Berger 
sphere to collapse at both ends. We are looking for a t = tsoi > such that 

{f(tsol),f{tsol), h{tsoi),h{tsoi),u{tsoi),u{tsoi)) = (0, -1, h, 0, u, 0) 

for some h > and some u > —2. We output initial conditions where 

SOLit) = fit) + (fit) + 1)2 + h\t) + u\t) 

satisfies SOL < 0.005 for some t > 0. The plot below indicates the output of the program. The 
parameters searched are < h < 3, —2 < u <2 and the stepsize is 0.005. 
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Figure 1 : Output of soliton hunting program on CP2 blown-up at 1 point 




Initial h 



The cluster about {h,u) = (0.7319,-0.5276) corresponds to the Koiso-Cao soHton, the cluster 
about (h, u) = (0.9595, 0) is the Page metric, and the third cluster is the Koiso-Cao soliton with 
the conjugate complex structure. 

We investigated for the soliton the behaviour of the winding angle of §1 numerically, and found 
that the angle decreases monotonically along the flow. 



5.2. Numerical results on Here G = 50(3) x S0{3) and the principal orbits are x S^. As 
we have mentioned before, Bohm [Bl| found infinitely many cohomogeneity one Einstein metrics 



on this manifold. We choose the normalisation | = —0.04 so that the standard round metric has 
initial conditions (/i(0),n(0)) = (10,0). 

We need one of the factors to collapse at one end, and the other factor to collapse at the other 
end. We are looking for a t = tsoi > such that 



ifitsol), f{tsol),Kisol), h{tsol),u{tsol),u{tsol)) = (/, 0, 0, -1, U, 0) 



for some / > and some u. We output initial conditions where 



SOLit) = Pit) + h^t) + (Ht) + 1)2 + ii^it) 



satisfies SOL < 0.005 for some t > 0. 



COHOMOGENEITY ONE SHRINKING RICCI SOLITONS; AN ANALYTIC AND NUMERICAL STUDY 21 



Figure 2: Output of Soliton Hunting Program on S5 




Initiai h 

The only values the program plotted are on the 'Einstein axis', u = 0. We can see a cluster 
of points around ^ = 10 which correspond to the standard metric, a cluster around 2.5 which is 
the first Bohm metric (h ^ 2.53554), then a final large cluster around (of course there should 
be infinitely many clusters between 2.5 and but obviously our numerics cannot detect all these). 
The authors are yet to find a cluster away from the Einstein axis. It is interesting to restrict the 
algorithm to = and output the value of SOL closest to for each h. This plot gives us 



Figure 3: Closest point on trajectory to closing up on S5 



O 
CO 




6 8 
Initial h 



We can clearly pick out the first three values of h corresponding to the Bohm metrics h 
10, 2.5354, 0.53054. We also see that the behaviour gets increasingly complicated as ^ — > 0. 
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5.3. Numerical results on x S^. As in the case, the group acting on the manifold is 
G = 50(3) X 50(3) and the principal orbits are x S^. 

We now need the same factor to collapse at each end. We are looking for a t = tgoi > such 
that 

{f{tsol),f{tsol), h{tsol),h{tsol), u{tsoi),u{tsoi)) = (0, -1, h, 0, u, 0) 

for h > and some u. As in the S"'' case, we do not find any non-trivial solitons. The hunting pro- 
gram does output clusters of points corresponding to the Bohm-Einstein metrics on this manifold. 
As it is more illuminating we output the smallest values of SOL along some trajectories with u = 0. 

Figure 4: Closest point on trajectory to closing up on 83x82 




1 2 3 4 5 6 

Initial h 



So we see we can pick out the standard product metric on this space h = 5 and then the Bohm 
metrics corresponding to the initial values h « 1.1779 and 0.23571. 

As mentioned earlier, the analytical methods used by Bohm to construct Einstein metrics only 
work in the dimension range 5 < dimM < 9. It is not clear, however, for which orbit types this 
restriction is essential. 

It is also conceivable that soliton solutions might exist on manifolds where Einstein metrics do 
not. This is certainly the situation in the Kahler case. 

We therefore examine next an example in higher dimensions, namely S^^. 

5.4. Numerical results on S-^^. In this example the group is G = SO{6) x 50(6) and the principal 
orbits are §^ x S^. 

The only cluster found by this search is one around (10, 0) corresponding to the standard round 
metric on S^^. 

5.5. Numerical results on HP^+^ft HP"^^- Another interesting case is that of the connected sum 
of two copies of quaternionic projective space with opposite orientations. In this case the groups 
are G = Sp{l) x Sp{n + 1), K = ASp{l) x Sp{n) and H = Sp{l) x Sp{l) x Sp{n). The principal 
orbits are G/K = S^""*"^, and we have di = 3,^2 = 4n and Cq = 4n + 8. The special orbits are 
MP", and the constant = 3. 
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In the case n = 1 Bohm [Bl] was able to prove the existence of a cohomogeneity one Einstein 

metric on MP^PP . He was also able to give numerical evidence for the existence of two coho- 
mogeneity one Einstein metrics in many cases where n > 2, although these examples are outside 
the dimension range where his analytical methods work (some of these examples were also found 
numerically by Page and Pope |PPj ). 

We also looked for solitons in these cases. We ran the hunting program for a variety of values of 
n, using the initial values h G (0.1, 10) and u £ (— ^^^^). In order to increase accuracy, we 
used a step size of 0.001 in the Runge-Kutta algorithm. 

We could not find any non-trivial solitons. When n = 1 we found the Bohm metric at /i ~ 
1.060793. We include a plot of SOL below for this case. It demonstrates a certain instability that 
all the Einstein metrics found on these manifolds shared, namely that one needs to be close to the 
exact conditions in order to get the trajectories to come close to closing up. 



Figure 5 : Closest point on trajectory to closing up on connect sum of HP2 




1 1.1 1.2 1.3 1.4 1.5 1.6 

Initial h 



We also include the same plot when n = 2. Notice we recover metrics that Bohm found numeri- 
cally with h ^ 1.0856062 and h ^ 0.2184791. 
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Figure 6: Closest point on trajectory to closing up on connect sum of HPS 




0.2 0.4 0.6 0.8 1 1.2 1.4 

Initial h 

5.6. Numerical results for F""*"^. We also considered the space 

pn+i £^Qj^ Yoi this space 

G = Sp{n + 1),K = U{1) X Sp{n) and H = Sp{l) x Sp{n), so the principal orbit is Cp2"+\ and 
the special orbits are MP". The fibration H/K ^ G/K ^ G/H is the twistor fibration for MP", 
and we have di = 2 and d2 = 4n, as well as Cq = 4n + 8 and = 8. 

Again we can recover Einstein metrics that Bohm and Gibbons-Page-Pope [GPPj found numer- 
ically, but we did not find any nontrivial solitons. We include below the results for n = 1 where we 
can see there is evidence for an Einstein metric at /i ~ 0.866. 

Figure 7: Closest point on trajectory to closing up on F2 




Initial h 

5.7. Numerical results on CaP^jjCaP . Finally we also considered the connected sum of two 
copies of the Cayley plane with opposite orientations. Now G = Spin{9),K = Spin{7) and H = 
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Spin{8), so the principal orbit is S^^, and the special orbits are We have di = 7 and d2 = 8, 
with RicQ = 28 and \\A\\^ = 7. 

We ran the soliton hunting program but did not find any evidence of Einstein metrics or solitons. 

5.8. A Non-compact example - the Gaussian. We conclude by discussing the smooth non 
compact Gaussian of Example 14.111 We consider the case when the principal orbits are x . 
With our normalisations that C = in equation (j5.4p and | = —4, the Gaussian soliton is described 
by 

dt^ + t^dSl + ]^dSl 

where is the usual Einstein metric on with Einstein constant 1. The soliton potential 
function is then given by u{t) = 2t^ — |. When we input the initial condition = 0.5 and u = —1.5 
and integrate we indeed recover the Gaussian soliton for a short time. However, as the numerical 
errors accumulate along the trajectory the numerical solution deviates from the Gaussian and very 
quickly becomes singular. Initial conditions close to the Gaussian also display this behaviour. This 
suggests instability of the soliton, so that it might be difficult to produce complete non-compact 
solutions close to the Gaussian. It should be noted that there are some gap theorems for the 
Gaussian soliton on Euclidean space, due to O. Munteanu-M.-T. Wang |MuW] and Yokota [Y] . 
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